J Arid Land 
https://doi.org/10.1007/s40333-018-0098-7 


Planting density affected biomass and grain yield of 
maize for seed production in an arid region of 
Northwest China 


JIANG Xuelian!’, TONG Ling'’, KANG Shaozhong', LI Fusheng*, LI Donghao!, QIN 
Yonghui', SHI Rongchao’, LI Jianbing* 


' Center for Agricultural Water Research in China, China Agricultural University, Beijing 100083, China; 

? Key Laboratory of Biochemistry and Molecular Biology in Universities of Shandong, Weifang University, Weifang 261061, 
China; 

> College of Agriculture, Guangxi University, Nanning 530005, China; 

4 Environmental Science and Engineering Program, University of Northern British Columbia, Prince George BC V2N 4Z9, 
Canada 


Abstract: Field experiments were conducted from 2012 to 2015 in an arid region of Northwest China to 
investigate the effects of planting density on plant growth, yield, and water use efficiency (WUE) of maize 
for seed production. Five planting densities of 6.75, 8.25, 9.75, 11.25 and 12.75 plants/m? were conducted 
in 2012, and a planting density of 14.25 plants/m? was added from 2013 to 2015. Through comparison with 
the AquaCrop yield model, a modified model was developed to estimate the biomass accumulation and yield 
under different planting densities using adjustment coefficient for normalized biomass water productivity 
and harvest index. It was found that the modified yield model had a better performance and could generate 
results with higher determination coefficient and lower error. The results indicated that higher planting 
density increased the leaf area index and biomass accumulation, but decreased the biomass accumulation 
per plant. The total yield increased rapidly as planting density increased to 11.25 plants/m?, but only a slight 
increase was observed when the density was greater than 11.25 plants/m?. The WUE also reached the 
maximum when planting density was 11.25 plants/m*, which was the recommended planting density of 
maize for seed production in Northwest China. 
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1 Introduction 


Planting density is one of the most important agronomic factors affecting crop grain yield and water 
use efficiency (WUE). An increased planting density usually improves the seasonal interception of 
solar radiation (Harper et al., 1979; Papadopoulos and Pararajasingham, 1997; Westgate et al., 
1997) that leads to increased canopy photosynthesis and biomass accumulation (Loomis and 
Connor, 1992; Coetto et al., 2013), and thus a higher grain yield and water productivity would be 
reached (Lang et al., 1956; Holt and Timmons, 1968; Fulton, 1970; Qiu et al., 2013). However, the 
grain yield tends to decline when planting density is above a certain level (Mohamed, 1999; Griesh 
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and Yakout, 2001) because of limited supply of carbon and nitrogen under intense interplant 
competition for intercepted radiation, soil nutrient and water. This would also increase the 
barrenness and decrease the kernel number (Lemcoff and Loomis, 1986). Consequently, an optimal 
planting density exists for maximizing the utilization of available resources and achieving the 
maximum yield per unit area by coordinating crop population and individual development. For 
example, a few studies have reported that maize grain yield and WUE increased with planting 
density when it was below an optimum level (Lang et al., 1956; Holt and Timmons, 1968; Westgate 
et al., 1997). The optimal planting density of maize varies depending on maize varieties and 
microclimate conditions. Olson and Sanders (1988) noted that the optimal maize planting density 
ranged from 3.00 to over 9.00 plants/m?. Under the surface irrigation conditions, the optimal 
planting density was found from 5.00 to 5.60 plants/m? in Egypt (Mohamed, 1999; Griesh and 
Yakout, 2001), while El-Hendawy et al. (2008) observed an optimal planting density of 7.10 
plants/m? under full irrigation. 

In recent years, the planting area of maize has increased rapidly in China, while the demand for 
maize seeds has also significantly increased. Particularly, in the arid region of Northwest China, 
the natural isolation conditions, abundant resources of light and heat are quite suitable for maize 
for seed production. As a result, the maize for seed production has become the main irrigated crop 
in this region with a rapidly increased planting area. As different from hybrid maize, the maize for 
seed production includes the female and male parents of short inbred lines without heterosis. A 
variety of studies have been conducted to examine the effects of planting density on the yield of 
common maize (Fulton, 1970; Mohamed, 1999; Griesh and Yakout, 2001; El-Hendawy et al., 2008). 
However, few studies were reported to investigate the effects of planting density on the yield and 
WUE of maize for seed production, especially by using yield models. The AquaCrop model 
developed recently by the Food and Agriculture Organization of the United Nations (FAO) assumed 
that the growth rate of aboveground biomass was linearly proportional to transpiration through 
biomass water productivity (WP), while the yield was determined by the harvest index (HI) and 
multiplication of the aboveground biomass (Raes et al., 2009; Steduto et al., 2009). It was found 
that HI increased as the increase of planting densities but started to decline after the optimal 
planting density achieved (DeLougherty and Crookston, 1979; Rahmati, 2009). To obtain a 
conservative estimation for a given crop, we should normalize WP for climate by taking into 
account the ratio of transpiration to reference crop evapotranspiration (Stedutoet al., 2009). For 
example, the normalized biomass water productivity (WP*) was found to be about 30 to 35 g/m? 
for C4 crops (Raes et al., 2011). However, it is unclear whether WP* and HI in the AquaCrop yield 
model can be taken as constants under different planting densities of maize for seed production, 
and this needs further investigation. The objectives of this study were then (1) to investigate the 
effects of planting density on plant growth, yield and WUE of maize for seed production; (2) to 
modify WP* and HI in the AquaCrop yield model and improve the estimated accuracy under 
different planting densities; and (3) to determine the recommended planting density of maize for 
seed production in the arid region of Northwest China. 


2 Materials and methods 


2.1 Experimental site 


The experiments were conducted during the growing seasons from 2012 to 2015 at Shiyanghe 
Experimental Station of China Agricultural University (37°52’N, 102°50'E; 1581 m a.s.l.) located 
in Wuwei, Gansu Province, China. This station is in a typical continental temperate climate zone, 
with the mean annual precipitation of 164.4 mm, mean annual pan evaporation of 2000 mm, 
average groundwater table lower than 25 m below the ground surface, mean annual sunshine 
duration over 3028 h, mean frost-free days of 150 d, and annual mean temperature of 8.8°C (Jiang 
et al., 2014). The soil texture at the experimental station was a light sandy loam with the mean soil 
dry bulk density of 1.38 g/cm? and field water capacity of 0.29 cm?/cm?. 


2.2 Experimental design 
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Maize for seed production is normally planted at a planting density of 9.75—-10.5 plants/m? in 
Northwest China. In 2012, the experiments were designed five planting densities, including 6.75, 
8.25, 9.75, 11.25 and 12.75 plants/m? referred as D1, D2, D3, D4 and D5, respectively. A planting 
density of 14.25 plants/m? was added from 2013 to 2015, referred as D6. A randomized complete 
block design with three replicates per planting density treatment was used. The experimental plot 
size was in 9.6 m long and 6.0 m wide. Maize in all plots was sown in each one row-male parents 
with seven-row female parents in 2012, and with five-row female parents from 2013 to 2015. Each 
planting density treatment had the same row spacing of 40 cm, and thus the six different planting 
densities of D1, D2, D3, D4, D5 and D6 were associated with different plant spacing within row of 
37, 30, 25, 22, 20 and 18 cm, respectively. Before sowing, all plots were fertilized with N of 136 
kg/hm?, P205 of 240 kg/hm? and KO of 50 kg/hm? as basal fertilizers. After fertilization, the soil 
surface in each plot was partly covered with 1.2 m width of film and 0.4 m width of bare soil 
between two rows. The plots were top-dressed with N of 364 kg/hm? on 10 June 2012, 5 June 2013, 
8 June 2014, and 8 June 2015, respectively. The irrigation quota was 100 mm using a border 
irrigation method. Table 1 lists the meteorological variables, irrigation and planting time of maize. 


Table 1 Meteorological variable as well as the planting and irrigation time over the whole growth stage of maize 
for seed production in 2012-2015 


Meteorological variable Planting and irrigation time 
Year R, È VPD Ta £ one 
7 Planting Irrigation 
(W/m?) (mm) (kPa) (°C) 

2012 265.38 129.40 1.23 19.02 16 Apr (Female) 6 Jun 
23 Apr (Male) 26 Jun 
26 Apr (Male) 13 Jul 
— 8 Aug 
— 27 Aug 
2013 208.73 68.40 1.45 19.29 13 Apr (Female) 6 Jun 
20 Apr (Male) 26 Jun 
23 Apr (Male) 13 Jul 
— 8 Aug 
— 27 Aug 
2014 222.06 206.2 1.28 17.94 15 Apr (Female) 1 Jun 
22 Apr (Male) 1 Jul 
25 Apr (Male) 20 Jul 
— 23 Aug 
2015 228.78 142.00 1.46 19.22 17 Apr (Female) 8 Jun 
24 Apr (Male) 2 Jul 
29 Apr (Male) 21 Jul 
— 7 Aug 
— 26 Aug 


Note: Rs, solar radiation; P, precipitation; VPD, vapor pressure deficit; Tą, air temperature; —, no data. 


2.3 Measurements 

2.3.1 Meteorological variables 

The meteorological variables including solar radiation (Rs), precipitation (P), air temperature (Ta) 
and relative humidity (RH) during the growing seasons of 2012-2015, were continuously observed 
by using a standard automatic weather station (Hobo, Onset Computer Corp., USA) that was 100 
m away from the experimental plot. The 30 min averages of all meteorological variables were 
calculated and recorded using a data logger. 

2.3.2 Leaf area index (LAI) and biomass accumulation 

Five female plants in each plot were randomly selected to measure the LAI and biomass 
accumulation. Leaf length and maximum width were measured with a ruler at intervals of 7—10 
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days after 15 days since sowing. The leaf area was determined by summing the rectangular area 
(lengthxmaximum width) of each completely developed leaf, and then adjusted by a factor of 0.74 
(Li et al., 2008). The leaf area per plant was divided by the surface area per plant to obtain LAI. 
The plant samples were collected at each growing stage. The total biomass accumulation was 
obtained after the leaves were oven dried at 60°C when the weights of plant samples kept constant 
in three hours, and then weighed using an electronic scale with the precision of 0.01 g. 


2.3.3 Determination of transpiration 


Soil water content in the root-zone was monitored by Diviner 2000 system (Sentek Pty Ltd., 
Australia). Two PVC access tubes were installed below the mulching and bare soil respectively in 
each plot. The measurements were made at 0.1 m intervals within 0—1 m soil depth every 5-7 days 
as well as before and after each irrigation and heavy rain event. The calibration of measurements 
was conducted using gravimetric method as near as possible to the probe. 

The evapotranspiration (ET (mm)) was estimated through soil water balance analysis (Rana and 
Katerji, 2000). Since the experimental plot was flat and the rainfall was not intensive, the surface 
runoff was neglected. The groundwater recharge was also negligible since groundwater level was 
lower than 25 m below the ground surface. Moreover, the measured soil water content at the 90- 
100 cm soil layer did not change before and after each irrigation event, thus the drainage can be 
ignored. ET is then estimated by Equation 1. 


ET=P.+I-AW, (1) 
where Pe is the effective precipitation (mm), I is the irrigation amount (mm), and AW represents 
the water content change (mm) in the root-zone that is determined by Equation 2. 


AW =W, -W,, (2) 


where W, and W, represent the mean water content (mm) in the root-zone at time f and tz, 
respectively. 

The soil evaporation (E (mm)) was measured with two micro-lysimeters placed within the bare 
soil between two plastic films in each plot. The micro-lysimeter with 20 cm height was made from 
PVC tubes with the inner and outer cylinder diameters of 10 and 11 cm, respectively. The outer 
cylinder was fixed into the soil with its top leveling with soil surface. The inner cylinder inside the 
outer cylinder was filled with an intact soil core, and was weighed at 19:00 LST every day by an 
electronic scale with the precision of 0.1 g. After the measurement of E, the crop transpiration under 
different planting densities during the growing season (T,(mm)) was calculated by Equation 3. 


T =ET-E. (3) 
2.3.4 Grain yield, water use efficiency (WUE) and harvest index (HI) 


Fifteen plants in the center of each plot were randomly selected for manual harvesting in each 
season. The seeds were weighed after sun-drying to obtain grain yield (Y, t/hm?). The WUE (kg/m?) 
was calculated using Equation 4. 


Y 
WUE = —. (4) 
ET 
The HI (%) under different planting densities was given by Equation 5. 
Y 
HI=—, (5) 
B; 


where Ba is the aboveground biomass accumulation (t/hm?). 
2.4 Estimation of biomass accumulation and yield 


In AquaCrop model, the Ba was estimated using the crop transpiration under different planting 
densities during the growing season (7;) under different planting densities during the growing 
season and normalized water productivity under standard crop management practice (WP, ) by 
Equation 6 (Steduto et al., 2009): 
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T 
ET, 


B, =K, WP; >, (6) 
where ETo is the reference crop evapotranspiration calculated according to the FAO Penman- 
Monteith equation (Allen et al., 1998); Ks, is the air temperature stress coefficient which can be 


taken as 1 due to no temperature stress occurred during the growing season in this study; WP, is 
the normalized water productivity under standard crop management practice which was represented 
by a reference planting density of 11.25 plants/m? in the study area. 

The harvestable yield (Ya (t/hm?) was then calculated by Ba and HI under standard crop 
management practice (HIo) using Equation 7. 


Y. = HB, (7) 
The AquaCrop yield model was modified in this study by introducing the adjusted WP* under 
different planting density (WP,,,) and the adjusted HI under different planting density (HIaaj). A 


linear relationship (Eq. 8) was assumed between the relative WP* (WP..,/W P} ) and the relative 
planting density (Daaj/Do): 


WP D. 
Aah, (8) 
WP, Dy 


where Daqj is the planting density (plants/m’) and Do is the reference planting density (11.25 
plants/m7). a and b represent the empirical coefficient which can be fitted using the measured 
transpiration and aboveground biomass accumulation in 2012 and 2013, and they were found as 
0.37 and 0.61, respectively. HIaaj can be calculated as Equation 9. 


-HI 


adj O 
=cD adj 


+d, (9) 
adj 0 
where c and d represent the empirical coefficients which were fitted using the measured yield and 
aboveground biomass accumulation in 2012 and 2013, and they were found as —0.18 and 1.29, 
respectively. 
The WPi; and Hlaaj can be calculated by Equations 10 and 11 based on Equations 8 and 9. 


* D * 
WP, -|2 = 10 WP (10) 


HI =HI, +(cDy,+d)(Dy, -D,). (11) 


The aboveground biomass accumulation (Ba) and grain yield (Ya) under different planting 
densities were then obtained using the modified model (Eqs. 12 and 13). 


Dag P T. 
Boa t we y À (12) 
¥,=(HI,+(cD,y+d)(D,y-D,)) B: - (13) 


2.5 Data analysis and evaluation of model performance 


SPSS 13.0 software (SPSS Inc., USA) was used for statistical analysis. The means were compared 
using Duncan's multiple-range test at the 5% probability level. The modeling performance was 
evaluated based on a linear regression between the estimated (E;) and observed (Qj) values of yield. 
Meanwhile, the mean absolute bias error (MAE) and root mean square error (RMSE) were included 
and were calculated as Equations 14 and 15 (Legates and McCabe, 1999), respectively. 


1 N 
MAE=— > — E, 
N >lo, ! 


; (14) 
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vs 2) ¥2 
Rost=| SCE -0,) | : (15) 
i=1 
where N is the number of data samples. 


3 Results and discussion 


3.1 Effect of planting density on LAI and biomass accumulation 


Figure 1 shows the LAI of maize for seed production under different planting densities in all 
seasons. Under each planting density, the LAI started to increase rapidly at the shooting stage, 
reached the maximum at the heading stage, and declined at the maturity stage due to leaf yellowing 
or wilting. During each season, the LAI increased with planting density, except at the seeding stage. 
The LAI difference between different planting densities reached the maximum at the heading stage. 
The maximum LAI differences of about 60% in 2012, 70% in 2013, 88% in 2014, and 80% in 2015 
were observed between D1 planting density and the maximum planting density, respectively. 
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Fig. 1 Leaf area index (LAI) of maize for seed production under different planting densities. Error bars denote 
standard error of the mean. D1, D2, D3, D4, D5 and D6 are planting densities of 6.75, 8.25, 9.75, 11.25, 12.75 and 
14.25 plants/m?, respectively. 


Figures 2 and 3 present the variation of biomass accumulation and the biomass per plant with 
planting densities in different years. It was observed that the biomass accumulation and biomass 
per plant showed a logistic growth pattern, and increased gradually with the advance of plant 
growth regardless of planting density. Except at the seeding stage, the biomass accumulation 
increased rapidly as planting density increased to 12.75 plants/m? (D5), followed by only slight 
increase with planting density (Fig. 2). However, higher planting density had an adverse effect on 
the biomass per plant in all seasons (Fig. 3). The biomass accumulation difference and biomass per 
plant difference between the lowest and highest planting densities reached the maximum at the 
maturity stage. For example, compared to the lowest planting density (D1), the highest planting 
density increased biomass accumulation by 55% in 2012, 67% in 2013, 60% in 2014 and 67% in 
2015, but decreased biomass per plant by 18% in 2012, 21% in 2013, 24% in 2014 and 21% in 
2015, respectively. This can be explained by the fact that the aboveground biomass accumulation 
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was associated with LAI which increased with planting density. A higher LAI increased the 
proportion of available radiation interception which then led to higher canopy photosynthesis and 
biomass accumulation (Loomis and Connor, 1992; Coetto et al., 2013). However, higher planting 
density resulted in more interplant competition, and thus led to reduced biomass per plant. The 
above results are consistent with previous findings (Akmal et al., 2014; Thimmappa et al., 2014; 
Dai et al., 2015). 
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Fig.2 Aboveground biomass accumulation of maize for seed production under different planting densities. Error 
bars denote standard error of the mean. 
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Fig. 3 Biomass per plant of maize for seed production under different planting densities. Error bars denote 
standard error of the mean. 
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3.2 Effect of planting density on yield and WUE 


Table 2 presents the yield and WUE results of maize for seed production under different plant 
densities from 2012 to 2015. It was found that the results showed a large difference in the study 
years due to the difference in microclimate conditions. The yield per plant and grain yield in 2012 
were higher than in other years, while the yield per plant declined with planting density in all the 
study years. Similar results were also found by Griesh and Yakout (2001) and El-Hendawy et al. 
(2008). This result might be attributed to more interplant competition at high planting density which 
led to reduced biomass per plant (Fig. 3) and thus a decreased yield per plant (Kamel et al., 1983; 
Soliman et al., 1995). 

However, the grain yield illustrated different variation trend from yield per plant. It increased 
rapidly as the planting density increased to 11.25 plants/m* (D4), followed by a slight increase 
when the planting density was above 11.25 plants/m*. The yield started to decrease when the 
planting density was above 12.75 plants/m?. This can be explained by the fact that under lower 
planting density (D1—D4), the interplant competition was small and thus the grain yield increased 
with the radiation interception by leaf area which increased during reproductive growth (Johnsonb 
et al., 1982; Duncan, 1986). However, when the planting density was great enough to intercept 
essentially all of the radiation at full canopy, further increase in planting density will not increase 
the yield (Duncan, 1986). Moreover, a higher planting density may lead to the formation of fewer 
flower initials, poor pollination resulting from asynchrony of tasseling and silking, or the abortion 
of kernels after fertilization, and consequently, a reduction of grain yield (Karlen and Camp, 1985; 
Hashemi-Dezfouli and Herbert, 1992). 

The WUE ranged from 1.84 to 2.46 kg/m? in 2012, 0.90 to 1.19 kg/m? in 2013, 1.13 to 1.39 
kg/m? in 2014 and 0.86 to 1.36 kg/m? in 2015, respectively (Table 2). The WUE in 2012 was higher 


Table 2 Yield and water use efficiency (WUE) of maize for seed production under different plant densities 


Yield per plant Grain yield ET WUE 
Year Treatment 
(g/plant) (t/hm?) (mm) (kg/m?) 
D1 121.2267 8.024 435.74° 1.844 
D2 116.76% 9.77° 453.434 215° 
2012 D3 111.98 11.09° 473.97° 2.34? 
D4 106.95 12.22" 495.68? 2.46% 
D5 100.904 12.66* SGT 2.46% 
D1 56.40? 3.814 421.80! 0.90° 
D2 54.20% 4.47¢ 443.30° 1.01° 
D3 53.10% SFS 462.804 1m2 
2013 
D4 50.50% 5.68% 486.10° 1i 
D5 47.50° 6.05° 507.10° ils 
D6 38.804 SS 531,90? 1.04° 
D1 Tia 4.814 424.08! i) 1B 
D2 67.32? 555° 452.82° 23E 
D3 64.84° 6.32” 473.724 11,338 
2014 
D4 62.404 7.02? 503.60° 1.394 
DS 55.997 7.08° 519.79° 1.36° 
D6 49.30! 7024 542.79° 1.29> 
D1 73.208 4.12° 477.00" 0.864 
D2 Tso" 4.944 485.35° 102° 
D3 69.64? 5.66° 489.534 iGS 
2015 
D4 68.18° 6.39° 495.29° 129° 
D5 63.96° 6.80* 501.12° 1.36% 
D6 47,534 5.64° 505.83" ill 


Note: Different lowercase letters indicate significant differences among different treatments for the same year at the P<0.05 level. ET, 
evapotranspiration; WUE, water use efficiency. 
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due to higher grain yield. However, the impact of planting density on WUE was similar for all the 
study periods. The WUE increased significantly with planting density when it was below 11.25 
plants/m? (D4), but no significant increase was observed when the planting density was beyond D4. 
Table 2 showed that significant decrease of grain yield and WUE was observed when the planting 
density was below 11.25 plants/m? or above 12.75 plants/m?. Thus the densities of 11.25 or 12.75 
plants/m? were optimal under the conditions of this study. Cox and Otis (1993) reported the 
maximum grain yield of maize at 7.41 plants/m*. Tetio-Kagho et al. (1998) noted that the grain 
yield of maize increased parabolically to the maximum yield of 1080 g/m? at about 10.00 plants/m?. 
Sangoi et al. (2002) reported an estimated optimum plant density of maize at about 8.50 plants/m?. 
El-Hendawy et al. (2008) indicated that the optimal density of maize in Egypt was 4.80 or 7.10 
plants/m?. The optimal planting density in our study was significantly higher than that of hybrid 
maize. The main reason was that the parent lines of maize for seed production were different from 
hybrid maize since the female and male parents were short inbred lines without heterosis, and the 
smaller interplant competition of parent lines led to more tolerance of higher planting density. When 
the planting density was 11.25 plants/m?, WUE reached the maximum. However, the WUE did not 
significantly increase when the planting density was too high. Moreover, a higher planting density 
resulted in large evapotranspiration and LAI (Table 2; Fig. 1), leading to the waste of water 
resources. It also becomes more difficult in efficiently carrying out agronomic practice under a 
higher planting density. Therefore, it is more economical to select the recommended planting 
density of 11.25 plants/m? for maize for seed production in Northwest China. 


3.3 The modification of yield model under different planting densities 


Figure 4 presents the estimated biomass accumulation and grain yield under different planting 
densities using the AquaCrop yield model (Ba and Ya) and the modified yield model. Compared to 
the measured biomass accumulation (B) and grain yield (Y), the results (Ba and Ya) from AquaCrop 
yield model were higher than the observed values when the planting density was below the 
reference planting density, but lower than the observation when the planting density was above the 
reference planting density (Fig. 4). The determination coefficient (R”), MAE and RMSE was 0.78, 
1.41 and 1.85 t/hm? for biomass accumulation, and 0.88, 0.62 and 0.50 t/hm? for grain yield, 
respectively (Table 3). The main reason is that the WP“ and HI for maize under different planting 
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Fig. 4 Comparison of modeling results with the measured value B and Y in 2012-2015. Ya, yield estimated by the AquaCrop yield 
model; Yq, yield estimated by the modified AquaCrop yield model; Ba, biomass accumulation estimated by the AquaCrop yield model; 
Ba, biomass accumulation estimated by the modified AquaCrop yield model. 
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densities were taken as constants in the AquaCrop yield model, which resulted in large error in 
model estimation. In fact, they were not constants. For example, WP* increased with planting 
density until the optimal density as shown in Figure 5. Moreover, HI increased with planting density 
when it was below the reference planting density (Fig. 5), but started to decline when the planting 
density was above the optimum planting density. Similar results were also found in other studies 
(DeLougherty and Crookston, 1979; Rahmati, 2009). In this study, the effect of planting density on 
WP” and HI was considered in the modified AquaCrop yield model. The estimated biomass 
accumulation and grain yield by the modified model were closer to the observed values (Fig. 4), 
with higher R?, lower MAE and RMSE as compared to the original AquaCrop model (Table 3). As 
a result, the modified AquaCrop yield model was considered to have a better performance in 
simulating biomass accumulation and grain yield under different planting densities of maize for 
seed production. 
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Fig. 5 Parameter variation of maize for seed production under different planting densities in 2012-2015. WP’, 
normalized biomass water productivity; HI, harvest index. 


Table3 Statistical analysis of modeling results from AquaCrop yield model (Baand Ya) and modified yield model 
(Baand Ya) in comparison with the measured values (B and Y) 


Item Regression equation R? MAE (t/hmô) RMSE (t/hm?) 
By=1.04B 0.78 1.41 1.85 
Biomass accumulation 
By=1.00B 0.98 0.47 0.56 
Ya=1.03Y 0.88 0.62 0.50 
Grain yield 
Ya=0.97Y 0.97 0.33 0.21 


Note: R?, determination coefficient; MAE, mean absolute bias error; RMSE, root mean square error. 


4 Conclusions 


A field experiment was conducted in 2012-2015 to examine the impacts of planting density on the 
biomass accumulation and yield of maize for seed production in Northwest China, and the 
AquaCrop yield model was modified by taking into account the effect of planting density on the 
WP’ and HI. The results indicated that increasing the planting density led to the increase of LAI 
and aboveground biomass accumulation, but significantly decreased the yield and biomass per 
plant. The grain yield and WUE increased rapidly as planting density increased to 11.25 plants/m?, 
but there was only a slight increase in the grain yield when the planting density was above 11.25 
plants/m?. It is thus more economical to select a recommended planting density of 11.25 plants/m? 
for maize for seed production in Northwest China. Within the range of planting densities used in 
the experiments, the original AquaCrop yield model produced a large error, with significantly 
greater estimation than the measured biomass accumulation and grain yield when the planting 
density was below the reference planting density. It had a lower estimation than the measured value 
when the planting density was above the reference planting density. However, the estimated results 
using the modified yield model were closer to the measured values, with higher determination 
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coefficient, and lower MAE and RMSE, indicating that the modified yield model by considering 
the effect of planting density had a better performance in estimating the biomass accumulation and 
grain yield of maize for seed production in the arid region of Northwest China. 
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